Domain-wall excitons and optical conductivity in one-dimensional Wigner lattices 
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Motivated by the recent finding that doped edge-sharing Cu-0 chain compounds like Na3Cu2 04 
and NagCusOio are realizations of one-dimensional Wigner crystals, we study the optical spectra 
of such systems. Charge excitations in ID Wigner crystals are described in terms of domain wall 
excitations with fractional charge. We investigate analytically and numerically the domain-wall 
excitations that dominate the optical absorption, and analyse the dispersion and the parameter 
range of exciton states characteristic for the long-ranged Coulomb attraction between domain walls. 
Here we focus on the Wigner lattice at quarter-filling relevant for Na3Cu2 04 and analyze in particular 
the role of second neighbor hopping t2 which is important in edge-sharing chain compounds. Large 
t2 drives an instability of the Wigner lattice via a soft domain-wall exciton towards a charge-density 
wave with a modulation period distinct from that of the Wigner lattice. Furthermore we calculate 
the temperature dependence of the DC-conductivity and show that it can be described by activated 
behavior combined with a dependent mobility. 

PACS numbers: 71.10.Fd, 78.20.Bh, 73.20.Qt, 71.35.-y 



I. INTRODUCTION 

At very low density an electron gas is expected to crys- 
tallize and to form a Wigner crystal, as in this limit the 
Coulomb interaction among electrons dominates their ki- 
netic energjsi. In the 70's Hubbard^ as well as Kondo 
and Yamaji^ suggested, expanding the considerations 
by Wigner, that the distribution of electrons in cer- 
tain tetracyanoquinodimethane (TCNQ) charge trans- 
fer salts may be controlled rather by the Coulomb in- 
teraction than by the kinetic energy {'^ band width), 
such that the electrons form a generalized Wigner lat- 
tice (WL) on the underlying TCNQ chain structure. 
This view suggests a strikingly different nature of charge 
excitations namely domain walls with fractional charge 
rather than particle-hole excitations as in usual metals 
and semiconductors^. This proposal, however, may be 
challenged on the grounds that the resulting periodic- 
ity of charge modulation can alternatively be explained 
by a Akp charge density wave (CDW)^*^ arising from 
an instability of the Fermi surface even in models with 
short-range interactions. In fact there is only a gradual 
crossover between the WL and the Akp CDW, thus there 
is no clear distinction between the two on the basis of the 
charge modulation period possible^. 

In a recent study of a new class of charge-ordered com- 
pounds Nai+a;Cu02^, which contain edge-sharing Cu- 
O chains, it has been suggested that the magnetic and 
thermodynamic properties of these compounds can only 
be explained in terms of WL formatiorJ. Edge-sharing 
chains consist of CUO4 squares just like the Cu-0 planes 
of high-Tc cuprates, but they are differently linked. The 
edge-sharing arrangement of CUO4 squares meets the WL 
criterion of small band width in an optimal way due 
to the almost 90° Cu-O-Cu bonds (Fig.l). Unexpected 
complexity is added because, apart from a small nearest- 
neighbor hopping matrix element ti , the second neighbor 
hopping t2 has to be considered which turns out larger 



as a consequence of the structure. While this unusual 
feature does not affect the classical WL order imposed 
by the Coulomb interaction, it changes the Fermi sur- 
face topology^*', and thereby allows to distinguish the 
WL from the CDW on the basis of the modulation pe- 
riod. The Nai+a:Cu02 compounds thus provide a first 
example where an unambiguous distinction between the 
generalized WL and a Fermi surface related ikp CDW is 
possible. 

The electron interaction driven Akp CDW has to be 
distinguished from the more familiar 2kF Peierls insta- 
bility which arises from a modulation of hopping matrix 
elements due to the coupling to periodic lattice distor- 
tions, i.e., leading to a CDW centered on bonds rather 
than on the ions^^. In contrast the Wigner lattice is 
based not on quantum but on classical energy consider- 
ations, namely which distribution of localized electrons 
has the lowest Coulomb interaction. It is quite remark- 
able though, that the periodicity of the WL coincides 
with that of the Akp CDW which emerges from a Fermi 
surface instability, i.e., from a pure quantum mechani- 
cal effect of strongly correlated electrons in models with 
nearest-neighbor hopping. 

Strictly speaking, at finite hopping t; (Z = 1, 2, • • • ) the 
WL is a quantum solid^^'^^. The kinetic energy causes 
virtual transitions to neighbor sites and leads thereby 
to a quantum mechanical smearing of the electron po- 
sitions. Hence electrons should be rather visualized as 
electron clouds that form a WL as consequence of the 
long-range Coulomb interaction. This delocalization of 
the electron in the WL is a pure quantum effect, i.e., 
controlled by the interplay of kinetic and interaction en- 
ergies. The quantum mechanical bluring of electron po- 
sition has considerable effects on the physical properties 
of Wigner lattices, as discussed in Ref. 9 in the context of 
edge-sharing chains in Na3Cu204 and NagCusOio com- 
pounds, where superexchange interactions and thereby 
the magnetic properties are strongly influenced by the 
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quantum nature of the WL. 

The aim of the present paper is to explore the charge- 
excitations of the ID WL, which are characterized as 
domain-wall excitations. We calculated the optical con- 
ductivity and its temperature dependence, as it is hoped 
that forthcoming experiments may provide further ev- 
idence for the WL nature of the electron structure of 
the Nai+a;Cu02 edge-sharing chain compounds. We use 
here both numerical, i.e., zero and finite- temperature 
diagonalization, and analytical methods to arrive at a 
deeper understanding of the nature of charge excitations 
in ID Wigner lattices at quarter-filling. In particular, 
we find that the long-range (repulsive) Coulomb interac- 
tion among electrons leads to exciton states below the 
domain-wall continuum, which appear as strong exciton 
absorption in the optical conductivity in the case of small 
hopping ti. Remarkably the second-neighbor hopping 
processes t2 are found to contribute strongly to the ex- 
citon dispersion. These processes lift the degeneracy of 
the exciton state and the lower branch gets soft at mo- 
mentum q = Tr/2 and leads to an exciton instability at 
a critical value t2,c- The CDW state beyond t2,c has a 
modulation period twice as large as that of the WL. The 
charge modulation in the CDW state is weak as com- 
pared to the WL state, as infcrcd from calculations of 
the static charge structure factor. 

Doped edge-sharing chains are also building blocks of 
the intensively studied system Sri4_2:Ca2:Cu2404i, the 
so called telephone number compound9i^4i^*i&. The com- 
posite structure of these materials consists of both ladder 
and chain structureaii. While originally the attention 
was directed toward the electronic properties of the lad- 
ders, because of their structural similarity to the high- 
Tc cuprates^®, more recently the number of papers re- 
porting information concerning the chains is increasing. 
The magnetic properties of these compounds which de- 
pend strongly on the doping are usually attributed to 
the chainsiSiSfliSi. Recently a quintupling of the chain 
unit cell in Sri4Cu2404i due to charge ordering below 
^ 200 K was reportedSSiS^. An additional complexity 
of these compounds is due to the exchange of holes be- 
tween chains and laddersi^, i.e., the doping concentration 
of edge-sharing chains is difficult to infer precisely. More- 
over an incommensurate modulation results from a misfit 
between the unit cells of ladders and chains. Work by van 
Smaalen24 and particularly a neutron scattering study of 
Braden at alm^ illuminate the subtle aspects of the in- 
terplay between modulations of chains and ladders. It 
has also been argued that the misfit between chains and 
ladders may modify the charge ordering and hence the 
spin structure on the chains^. Particularly remarkable 
is a study by Isobe et alm^, who succeeded in resolving 
the internal charge structure of the charge modulation of 
the compound [{Sro.o29Cao.97i)2Cu203]54[Cu 02)77 with 
a chain unit cell containing 77 Cu sites. Inspection shows 
that the charge pattern found in the structure analy- 
sis compares favorably with that expected for a general- 
ized Wigner lattice. Recently also the modulation of the 



charge density in the ladders of Sri4Cu2404i has been 
reported by Abbamonte et alm^. 

Wigner crystals are prime examples for strongly corre- 
lated states in the sense that electrons do the utmost to 
avoid each other in real space^- . In general strong corre- 
lations (i.e., large on-site interaction U) and the associ- 
ated reduction of kinetic energy are favorable for charge 
localization. As a consequence the long-range Coulomb 
interaction may become relevant in strongly correlated 
systems, i.e., leading to charge ordered states and WL 
order in higher dimensions at particular fillings. Exam- 
ples are the manganites at quarter-filling which reveal 
checker-board charge orderi^Si^ and the layered molecu- 
lar crystals of the BEDT-TTF typo32,33i34,35^ ^jgQ charge 
stripes in high-Tc compounds at 1/8 doping^S reflect the 
interplay of strong correlations and long-range Coulomb 
interactions. 

The outline of the paper is as follows: In Section II 
we describe the Hubbard- Wigner model for edge-sharing 
chain systems and introduce the corresponding spinless 
fermion Hamiltonian. Furthermore we analyse the re- 
sulting domain-wall interactions for both Coulomb and 
truncated interactions. In Section III we present diago- 
nalization results for the charge-excitation spectrum and 
analyse the emerging excitonic states in the case of the 
model with Coulomb interaction and nearest-neighbor 
hopping. Here we also provide an analytical derivation 
both for the continuum and the exciton states. The tem- 
perature dependence of structure factor, kinetic energy, 
and optical conductivity calculated by means of exact 
diagonalization (ED). Section IV deals with the signifi- 
cant changes of the excitonic states introduced by second 
neighbor hopping. These are analysed with help of the 
analytical solution. In particular we show here that the 
soft exciton states near q = tt/2 are visible in the optical 
conductivity at elevated temperature as mid-gap absorp- 
tion. Finally the DC-conductivity of ID Wigner lattices 
is discussed in Section V, while our conclusions are sum- 
marized in Section VI. 



II. MODEL 

A. Wigner lattices in doped edge-sharing chain 
compounds 

As in the high-Tc cuprates Cu^"'' is in a (P configu- 
ration with spin 1/2, while Cu'^"'" is in a d^-ligand hole 
{cfiLh) singlet state, i.e. involving the O neighbors of 
the Cu ion, well known as Zhang-Rice singlet^. Yet in 
contrast to the high-Tc superconductors the edge-sharing 
geometry with nearly 90° Cu-O-Cu bonds (Fig^l leads 
to strongly reduced hopping matrix elements. This sets 
the stage for the long-range Coulomb force as predomi- 
nant interaction 

Hcoui = U ^7ii,|7ii,x + ^ ViniUi+i, (1) 

i i,l>l 



3 



Op,: 



Pv-7 



FIG. 1: Orbital structure of edge-sharing copper-oxygen 
chains. The 90 degree Cu-O-Cu hopping ti and the sec- 
ond neighbor Cu-O-O-Cu hopping path t2 are marked by ar- 
rows. Shading indicates the p-d covalent mixing, as well as 
the hole distribution in the d^Lh Zhang-Rice singlet states in 
Na3Cu204 with alternating charge order (dark shading). 



where the on-site interaction U takes care of the strongly 
correlated character of these systems and suppresses 
charge fluctuations leading to Cu^+ {d^^) configurations. 
In our model we associate the dPLh {(f,d^'^) ioniza- 
tion state with (1,2) electrons, respectively, and Ui^cr 
{a =t,i) counts the number of electrons with spin a, 
while rii — Thus the concentration of electrons 

p is related to the concentration of holes 6 — 1—p (relative 
to the configuration) used in the high-Tc literature^S. 
The Coulomb interaction Vi is screened by the polariza- 
tion of neighbouring chains and by core electrons^. We 
shall not try to explore the subtleties of screening due 
to the embedding, and assume for sake of simplicity a 
generic Coulomb law V/l = y, I = 1,2,...., and keep the 
nearest neighbor interaction y as a parameteraSE. Crucial 
for the following is that the interaction is long ranged and 
convex, i.e., V/' = Vi_i - 2Vi + Vi+i > 0. 

For commensurate doping concentration p = m/n 
the interaction Vi selects a particular charge ordering 
pattern^. The resulting charge order is immediately obvi- 
ous for the filling fractions p = 1/4, 1/3 and 1/2 (Fig|2Ia- 
c)), which involve an equidistant arrangement of the 
Cu^"*" sites (arrows in Fig|2I). For a general ratio p — m/n 
this leads to complex structures with unit cell size n. 
In case of p = 2/5 and 3/5 we encounter in Fig|2d,e) 
the charge order observed for Sri4_j;Caa;Cu2404i and 
NasCusOiG^. Charge localization, however, is not per- 
fect in Wigner insulators as electrons still undergo virtual 
transitions to neighboring sites (Figl^ff)) in order to re- 
tain partially their kinetic energy . The energy of the 
lowest excitations and the impact of kinetic energy de- 
pend strongly on p = m/n. For example, the energy of 
the lowest excitation relative to the ground state Figl^c) 
is ~ V2 while the excitation for p = 3/5 shown in Fig. 
2(f) is ~ I5", about an order of magnitude smaller. Hence 
quantum charge fluctuations are more important in the 
latter caseS. 

To investigate the role of kinetic energy we explore the 
dynamics of electrons starting from the one-dimensional 
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FIG. 2: Wigner charge order resulting from Coulomb repul- 
sion and associated modulated Heisenberg spin structure for 
p = 1/4, 1/3, 1/2, 2/5 and 3/5 doping (a-e). The spin-1/2 of 
Cu^^ (arrows) is responsible for magnetism, whereas Cu'^"'' 
(red circle) is nonmagnetic. The spin arrangement is that 
expected for ferromagnetic Ji and antiferromagnetic J2 ex- 
change interaction. The charge unit cells (shaded) contain 4, 
3, 2 and 5 sites, respectively. The structures (c) and (d) are 
realized in Na3Cu204 and NagCusOio, respectively. A fluc- 
tuation of spin-position due to a low-energy charge excitation 
is shown in (f) for p — 3/5. 



Hubbard- Wigner model Hhw — Hcoui + HKm^ , where 

Hxin = - XI ^1^4+1, ^(^i,<y + (2) 



describes the hopping of an electron with spin a from 
site i to site i -\- I and vice versa; and ni^^ = c\^c^^. 
Due to the almost 90 degree Cu-O-Cu angle the hopping 
ti between nearest neighbor Cu sites results mainly from 
direct d—d exchange, while t2 originates from hopping via 
a Cu-O-O-Cu path^ (Fig. 1), leading to the remarkable 
fact \t2\ > \ti\. We adopt here as typical values ii ~ 63 
meV, t2 94 meV, derived from ab-initio band structure 
calculations for the Cu^+ edge-sharing reference system 
Li2Cu0^i. The hopping integrals are much smaller than 
our estimates for [/ ^ 3.8 eV and V ^ 1.5 eV. Thus these 
parameters suggest that the edge-sharing chains are well 
inside the WL regime. 

In particular 1^2 1 > l^il implies that the second- 
neighbor (antiferromagnetic) exchange integral J2 is 
large compared to the nearest-neighbor interaction Ji 
and the interchain couplings. In fact this is consis- 
tent with the magnetic properties of the Na3Cu204 and 
Na8Cu50io compounds^ and also with recent neutron 
diffraction data of the compound NaCu202 which con- 
sists of p = 1 edge-sharing chains^^. Magnetic excitation 
spectra determined by Raman spectroscopy*'^ and a high- 
field NMR studji^ on NaCu202 single crystals confirm 
these conclusions. 

It is evident that the spatial variation of charges in 
the WL at the same time implies a complementary 
arrangement for the spins and thus leads to spatially 
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modulated Heisenberg spin chains with varying distances 
among the spinsSii^. This new category of spin mod- 
els shows very different magnetic properties for differ- 
ent commensurabiHties p — m/n, as actually observed 
in the Sri4_a;Caa;Cu2404i compound3^2*2i and in the 
Nai+a;Cu02 systems. For example, Sri4Cu2404i has the 
commensurability p — 2/5 corresponding to the structure 
shown in Figl^Jd) and its ground state is determined by 
singlet pairsSSiSiiS^l. Thus its magnetic properties are 
quite distinct from those of Na3Cu204 {p = 1/2) and 
NagCusOio (p = 3/5)a. 

B. Spinless fermions and domain- wall interaction 

Our study of the charge structure and dynamics will 
be based on the spinless fermion version of the Hubbard- 
Wigner Hamiltonian H = Ht + He ■ 

ij>l ij>l 

The underlying assumption of spin-charge separation, 
i.e., the neglect of the effects of spin-degrees of freedom 
on charge correlations and excitations, can be rational- 
ized when starting from the ID Hubbard model which is 
Bethe ansatz solublc^^. That model has marginal spin- 
charge coupling^Si^Si^iiS and its charge excitations are 
described by free spinless fermionsS. The addition of the 
long-range Coulomb interaction leads then to the model 
of interacting spinless fermions^'*. We note, however, that 
in the ID WL there are additional mechanisms due to 
the long-range Coulomb interaction which influence the 
exchange interactions as discussed in Ref. 0, and which 
induce some coupling between charge and spin. Such ef- 
fects will be neglected here. 

In the following we shall focus on the p = 1/2 case. 
The ground state has perfect alternating charge order as 
shown in Fig.OfaO) for ti = 0. At finite hopping ti this 
ground state remains stable, yet domain wall pairs get 
mixed in due to quantum fluctuations. An elementary 
ti hopping process involves the interchange of (xO)— >(0x) 
pairs in Fig.O The charge excitations in Wigner lattices, 
caused e.g. by optical excitations, involve the creation of 
domain wall (DW) pairs. These DW's can move as a 
consequence of -processes of the kinetic energy. The 
role of t2 processes is distinct in the p — 1/2 case; they 
are blocked in the perfect WL and contribute only in the 
presence of DW pairs. We shall consider this problem in 
a later Section, and assume for the moment ^2 ~ 0. 

The charge of a DW can be invoked from a Gedanken 
experiment due to Hubbard^. The excitation generated 
by adding an extra electron with charge e as shown 
in Fig. ISIbO) will dissociate into two equivalent DW's 
with fractional charge e/2 separated by regions of perfect 
charge order. In general the domain wall charge depends 
on the commensurability p ~ m/n. 

The creation of a domain wall pair requires an energy 
cx Vi. Domain walls can propagate freely, yet due to the 
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FIG. 3: (a) Charge excitations in a Wigner lattice move in 
form of domain walls (DW's). Here we consider the p = 0.5 
case where Cu^"'' ( Cu^"*") are indicated by x(o), respectively, 
(b) The addition of an extra electron with charge e leads to a 
high energy state which decays into two DW's with fractional 
charge e/2. 

long-range interaction Vi they attract each other. For 
p = 1/2 the energy of two domain walls at distance 
d ~ 2m is determined by the recursion relation 

OG 

A™ = A„_i +Y,VZ, m = 2,3, • • • (4) 

n— m 

with Al = y^'n- Here V/' - I^„i - 2Vi + Vi+i de- 

notes the discrete second derivative of the interaction. 
In the following we shall assume Coulomb interaction 
Vi — V/l, and furtheron use V = 1 as unit of energy. 
In this case Ai = 21n2 — 1. 

An excellent asymptotic expansion for the domain-wall 
interaction (DWI) has been derived by Fratini et al.^^ for 
the Coulomb case: 

A,„ « 1/2- l/8TO-f l/(4m)3 . (5) 

Remarkably this expression has an accuracy of two digits 
even at m = 1. It is interesting to note, that the leading 
interaction term — l/8m in the asymptotic expansion can 
be interpreted in terms of an effective Coulomb interac- 
tion between the fractional charges of DW's qi^2 = ±1/2. 
That is, the DW interaction is given as qiq2/d = —l/8m, 
where d = 2m is the distance between the DW centers. 
Hence the interaction between the domain walls A^ pro- 
vides a manifestation of the fractional charge of the do- 
main walls. 

It is evident that a truncated Coulomb interaction, i.e. 

= for Z > Imax, may not stabilize the Wigner lattice 
structure at general rational fillingsSi^SiS. The required 
value for Imax increases with the commensurability n, 
where p — m/n is the filling fraction. In the case of 
the excitations truncation is even worse. This problem 
is demonstrated in Fig. ^ which displays the energy A^ 
of two domain walls as function of distance d = 2m for 
the most simple Wigner lattice p = 1/2. In the case of 
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FIG. 4: Interaction energy Am of two DW's at distance 
d = 2m for a Wigner lattice with density p — 0.5. Result for 
Coulomb interaction (solid line) and asymptotic expression 
(circles) is compared with Am obtained for Coulomb interac- 
tion truncated at Imax ~ 1,2, ,6 (triangles and squares). 



Coulomb interaction A„j is an increasing function of m, 
i.e. the domain wall interaction is attractive, and con- 
verges against Aqo = 1/2. Remarkably truncation leads 
to two different classes of behavior, repulsion at short 
distance for Imax = even and attraction for Imax = odd. 
Note that in the former case A/ ^ for large I, while 
for Imax — 1 which is sufficient to stabilize the Wigner 
lattice for p = 1/2 there is neither attraction nor repul- 
sion, i.e., DW-pairs are not confined in this case. Hence 
calculations of excitation spectra based on models with 
truncated interaction Vi, V2, '^i^a^ O'" models with ar- 
bitraly chosen parameters Vi, V2 etc. must be considered 
with care when they are compared to the Coulomb case. 
A frequently studied model is the model with both near- 
est neighbor-hopping ti and nearest neighbor interaction 
Vi^i^SiS. The effects of longer-range interactions have 
been studied by Poilblanc et alm^, though for relatively 
weak interactions, i.e., outside the WL regime. 

In this work we shall confine ourselves mainly to the 
discussion of the long-range Coulomb interaction. For the 
analytical considerations we shall occationally consider 
the truncated models with Imax = 1 and Imax = 3. 



III. RESULTS FOR NEAREST-NEIGHBOR 
HOPPING 

A. Domain wall continuum and exciton 

As shown above charge excitations in Wigner lattices 
form domain walls (DW's) separating regions of perfect 
charge order. DW's move as a consequence of the kinetic 
energy operator Ht . For sake of transparency we consider 
first the model with nearest-neighbor hopping motion ti. 
Numerical results for the excitation spectra for Coulomb 
interaction are given in Fig. |S1 for a quarter-filled ring 



FIG. 5: Excitation spectrum for the quarter-filled chain 
with A' = 26 sites and Coulomb interaction (V — 1) for (a) 
ti/V = 0.02 and (b) ti/V = 0.05 as obtained by exact di- 
agonalization. The continuum due to domain-wall pairs is 
centered at iJ ~ 0.45y and its width is ~ 8ti. At small ti the 
bound state is well separated from the bottom of the contin- 
uum, while for larger ti (b) it merges with the continuum, but 
persists near q — tt/2. Parts of a continuum due to excitation 
of 4 domain walls is visible at the top of the figure. 



(i.e. p — 1/2) with A^ = 26 sites. The figure shows 
the two degenerate ground-states at 5 = and tt, the 2 
domain-wall continuum and part of the 4 DW continuum 
at high energy. The 2 DW continuum is centered near 
E ^ QA5V as expected from the DW interaction A„i « 

1/2 - l/8m + l/(4m)3 in the Coulomb case. The 

small downward shift from E — Aoo = 1/2 is attributed 
to the interaction with the 4 DW states. As a result of the 
Coulomb attraction between DW pairs an excitonic state 
emerges below the 2 DW continuum, which is expected 
to play a prominent role in the optical absorption. With 
increasing hopping ti the width of the 2 DW spectrum 
^ 8^1 increases and the exciton binding energy decreases, 
and eventually the exciton is absorbed by the continuum 
near q = and tt. Yet due to the dispersion of both the 
lower edge of the 2 DW continuum and the exciton, the 
excitonic state survives in the vicinity of q = 7r/2. 

Next we explore the analytical structure of the DW 
continuum and then analyse the dispersion of the exciton 
in more detail. We will address three different problems: 
(i) The DW continuum in the case that all A; are equal. 
This case pertains to the model where only the n.n. inter- 
action Vi is kept. There is no exciton state in this case 
and DW's are not confined^'*'^^. (ii) Next we explore 
the solution for the bound state emerging for Ai < A^o 
and A; — Aqo for I > 2. This case directly applies to 
the model where the Coulomb interaction is truncated at 
Imax — 3, but this solution also provides an approximate 
description for the case with full Coulomb interaction, if 
one adopts appropriate values for Ai and Aqq. Finally 
(iii), we analyse the special role of the second neighbor 
hopping t2 in the case p= 1/2. 
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We begin with the motion of a DW-pair as indicated 
in Fig. |3fa) and denote the pair state by |n,m >. Here 
n — m and n + m denote the centers of the xx and 00 
DW's indicated by vertical bars in Fig. O while n de- 
notes the center of mass coordinate of a DW pair. We 
shall consider periodic boundary conditions, that is, even 
numbered rings of size TV. Then it is useful to introduce 
the auxilliary Bloch states 



1 



n, rn); 



(6) 



where for periodic boundary conditions momenta are 

defined as — 2tiv/N and v — 0,1, , iV — 1. 

When applying the translational operator T/ one obtains 
Ti\il^q.m) = e~"''|'0</,m)- Thc DW interaction energy as- 
sociated with these states is A„, as discussed above: 



(7) 



The action of thc kinetic energy operator Ht on thc local 
DW states yields: 

Ht\n,m)^ti (\n-l , m + 1) (8) 

+ \n + l , m+ 1)^(1 -(5,„+i,jv/2) 

+ (\n+l , m- 1) 

+ \n-l , m - 1) )(1 - (5„i_i^o) • 



This can be expressed in terms of the auxilliary Bloch 
basis Eq.© as: 

Ht\lpq,„i) = ti{q) (1 - 5,„+l^Ar/2)|V'<?,rn+l) (9) 
+ (1 - Sm.-l,o)Hq,m~l) 

with ti{q) — 2ti cos{q). (i) In the case where all A„j have 
the same value, which we denote Aqo, the solution for the 
domain wall continuum is straightforward: 



N/2-1 



sin{pm)\'ipq 



(10) 



m— 1 



where the N/2 — 1 pscudo momenta are determined by 
Pf, = 27r/i/iV and fi = 1,2, ,N/2 - 1. The corre- 
sponding energies of the DW continuum are: 



+ 2ti[cos{q + p) + cos{q-p)]]. (11) 



That is, they are given as linear combinations of single 
domain wall energies. The total width of the continuum 
is 8ti. 

While the ti-Vi model with nearest-neighbor inter- 
action stabilizes the alternating charge-ordered ground 
state for p = 0.5, it does not lead to an attractive inter- 
action between domain walls. The excitation spectrum 
consists of the domain wall continuum. The absence of 
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Analytical solution for the domain wall contin- 
q^p of the quarter-filled chain with N — 26 sites and 



FIG. 6 

uum E, 

A„ = 1/2 for ti = 0.02 and 0.05. Exciton dispersion Eex{q) 
calculated for the model with Ai = 21n2 — 1 and A^ = 1/2 
for m > 2 relevant for the case of long-range interaction. 



the domain-wall exciton shows that domain walls are not 
confined in this caseSA*^. 

The simplest case which shows a bound state is the 
model where the Coulomb interaction is truncated at 
hnax — 3. The resulting DWI A^ in this case is Ai = | V 
and A,„ = Aqo = V for m > 2, i.e., the the DW con- 
tinuum is centered around V. Next we shall derive an 
analytical expression for thc exciton energy. 

The exciton solution can be obtained by expansion of 
the secular determinant 



D 



LUi ti{q) 
ti{q) UJ2 ti{q) 
ti{q) UJ2 ti{q) 



ti{q) LU2 h{q) 
ti{q) uji 



(12) 



where uji = Ai — E and ti(q) = 2ticos{q). 

D = UjjDM - 2uJih{qfDM-l + h{qfDM-2 (13) 



Dm = 



L02 ti{q) ... 
ti{q) UJ2 hiq) 



ti{q) UJ2 



Dm = UJ2DM-1 - ti{qYDM-2, 



(14) 



(15) 



where Di = UJ2, D2 — u>\ — ti{q)^ and _Do = 1. The last 
relation is required so that the previous determinantal 
equation is fulfilled for M > 2. The equation may be 
solved by the polynomial ansatz D]\i = p*^, leading to 



P 



1 r 

2 



UJ2 ± 



4h{qy 



(16) 



The equation for the bound state uip - 
the form: 



ti{q)'^ takes then 



( Ai -E)[A2-E± v/(A2~£;)2~4ti(g)2j = At, (9)' 

The physical solution, which is twofold degenerate be- 
cause of Ea. ((T^ . is 



Ai' 



(18) 



where the dispersion t{q) — 2tiCos{q), which appears 
squared, reflects the hopping of domain walls by two lat- 
tice units. We also notice that the dispersion is weighted 
by the DW binding potential A2 — Ai . As this quantity is 
y/3 in the model with truncated interaction {Imax = 3), 
the dispersion of the exciton (not shown) is smaller than 
in the model with Coulomb interaction. 

It is instructive to compare the analytical results 
obtained so far with the numerical data obtained for 
Coulomb interaction displayed in Fig. We approxi- 
mate A„j as follows Ai — 2 In 2 — 1 and A^ = 1/2 for 
TO > 2. The results are shown in Fig. El Although the ap- 
proximation used for the DWI Am is rather crude, there 
is quite a good agreement with exact diagonalization. 



B. Structure factor and melting of Wigner Lattice 

Wigner lattice order gets weakened with increasing ki- 
netic energy ~ ti. Yet in the model with only near- 
est neighbor hopping there is a gradual crossover to the 
CDW with weaker charge order but the same modula- 
tion periodSi, which is(7 = 7rat/5=l/2. This is the 2k p 
CDW of spinless fermions, and corresponds to the Akp 
CDW'^-S if spin is included in the model. The crossover 
from the WL to the CDW regime is reflected in the vari- 
ation of the charge gaps a,t q = tt{0) displayed in Fig[7Ka) 
as function of ti. While the two gaps are almost equal in 
the WL regime at small ti, for larger <i the q — tt{= 2kp) 
charge gap controls the CDW state. In this figure one can 
also see that in the WL regime the charge gaps calculated 
numerically for a = 26 site cluster are well described 
by the perturbative expression, Ea. H18|l . for the exciton 
energy E^.^ using Ai = 2 In 2 — 1 and A2 = 1/2 appropri- 
ate for the WL stabilized by Coulomb interaction. We 
note, that for the model with nearest neighbor interac- 
tion V the CDW solution for p = 1/2 is expected to be 
stable up to ti — 0.5V"—, at even higher values for ti a 
metallic state will be realizedSS. 

Next we shall investigate the disappearence of WL or- 
der due to thermal fluctuations. It is well known that 
one-dimensional systems controlled by short-range inter- 
actions do not exhibit long-range order at finite temper- 
ature. The corresponding absence of a phase transition 
at finite temperature follows from the Mermin- Wagner 
theoremSSiS. However, there is a finite transition tem- 
perature even in ID models if the decay of interactions 
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FIG. 7: (a) The variation of the charge gap at g = tt as 
function of t\ displays the crossover from the Wigner lattice 
to the CDW regime in the case with nearest-neighbor hopping 
t\. Results are obtained for a A'^ = 26 site ring using an 
energy cutoff 1.45 V . Exciton energy of WL (dashed curve) 
according to Ea. (ll8l l. (b) Dependence of the gap energies for 
large t\ — 0.1(0.25) on the cutoff energy. 



is power-law like and sufficiently slow, as has been shown 
for Heisenberg and Ising systems with ferromagnetio^ 
and antiferromagnetioS^ interactions. The latter work by 
Erwin and Hellberg which deals with an interaction de- 
caying as 1/R is of direct relevance for the WL case. 

The melting of the generalized WI— due to domain 
wall excitations is refiected in the structure factor N{q), 
whose temperature dependence we study here. The T = 
charge structure factor 



Niq) - (I/N! 



i^exp''(*-^)(n,nj 



(19) 



is normalized here with respect to the total electron num- 
ber of electrons A^e, such that N{q — n) — 1 for perfect 
CO, that is, at low temperature and for small ti. The al- 
ternating charge order is reflected in the peak of N{q) at 
q = TT which starts to decrease strongly for T > 0.051^. 
Yet even at much higher temperatures N(q) reveals a 
maximum at q — n which indicates the persistence of pro- 
nounced short range charge correlations. The results for 
N{tt) displayed in the inset in Fig.|S|indicate the melting 
of Wigner charge order^i at a temperature ~ 0.06^ , 
i.e., for the case of Coulomb interaction and ti = 0.021^. 

We note that T,n ^ 0.06V {V ^ 1 eV) has the correct 
order of magnitude as the charge order transition ~ 
455 K in the compound Na3Cu20^. Nevertheless we are 
far from any quantitative comparison, as 3D effects could 
lead to a further enhancement of or to a reduction due 
to frustration effects resulting, e.g., from the interaction 
with Na-ion potentials. 

Before we move on to the calculation of the optical con- 
ductivity, we briefly comment on the exact diagonaliza- 
tion scheme used. In the WL regime the eigenstates are 
energetically ordered according to the number of domain- 
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FIG. 8: Structure factor N{q) versus q for a A*' = 26 
chain with Coulomb interaction and ti — 0.02V for differ- 
ent temperatures. Inset shows the temperature dependence 
of iV(g = tt). 




FIG. 9: Optical conductivity poo"(<^) at T = (solid line) 
and density of states -0(0;) of 2- and 4-DW excitations (dashed 
line) for a A'^ = 26 chain with Coulomb interaction for ti/V = 
0.02 (top) and 0.05 (bottom), respectively. The spectra are 
shown with a Lorentzian broadening 7 = 0.03, while for the 
bound state in o-{uj) 7 — 10~* was used. The strong increase 
of D{uj) at a; ~ 0.7 marks the onset of the 4-DW continuum 
which does not contribute to 17(0;). 



wall pairs that are excited (c.f. FiglS)). This property of 
the Coulomb interaction, i.e., that it removes the high 
degeneracy of the usual Hubbard model, allows to trun- 
cate the Hilbert space by selecting configurations with a 
small number of domain walls, and to perform a full di- 
agonalization in the truncated space defined by a cutoff 
energy. This procedure is expected to work well in the 
WL regime, but it should become less accurate in the 
CDW regime. As the CDW can be visualized rather as 
a Fermi sea with a small gap due to the relevant 2kF 
scattering processes, a truncated domain wall basis will 
certainly not be appropriate or require an extremly large 
cutoff energy such that again the full Hilbert space is 
covered. Figure [TJb) shows the convergency of the low- 
est excitation energies at g = and tt as function of the 
cutoff energy. As can be seen from the figure, this trunca- 
tion works extremely well even for ti = 0.11^, i.e., at the 
crossover to the CDW regime. Thus most of our calcu- 
lations for the bigger iV = 26 clusters employ an energy 
cutoff 1.45. Also a basis of momentum eigenstates is used 
to further reduce the size of the matrices. 



C. Optical conductivity 

Optical conductivity experiments could provide impor- 
tant information about the size of the charge gap, the 
presence or absence of exciton features, and the shape 
and width of the domain-wall continuum. Thus in com- 
bination with the theoretical spectra optical data may 
allow to determine the basic parameters of the Hubbard- 
Wigner model more precisely. As we are dealing with 
insulating systems it suffices to focus on the finite fre- 
quency absorption (7{lu) as given by the Kubo equation 



in terms of the current-current correlation functionSS 



^-Ett^tt^^tfMttT' (20) 



Nlu 



where {n\ ((m|) are eigenstates with energy £'„ {E„i), re- 
spectively. The current operator jx for the lattice model 
is defined in the usual way as 



E 

i,i>i 



(21) 



where 61 = la is the hopping length, and a the Cu-Cu 
distance along the chain. The conductivity will be given 
in dimensionless formS^, po * with pQ = Hv/e^a^ 

where v — abc is the cell volume per Cu-site, and h/e^ = 
A.lkfl the von Klitzing constant. 

In this section we consider nearest-neighbor hopping 
amplitudes only. Equation |(5DI) can be evaluated in a 
straightforward way via exact diagonalization. Results 
for <j{uj) at T = are given in Fig. El for two different 
hopping matrix elements ti/V — 0.02 and 0.05, which re- 
veal the broadening of the domain-wall continuum ~ 8ti 
and the disappearance of the exciton peak at larger ti 
in the continuum. At ti = 0.02^ the exciton at energy 
0.38V has large weight and the shape of the DW con- 
tinuum is asymmetric and peaked at the lower frequency 
edge. The strong asymmetry of the continuum persists 
at larger ti, i.e., when the exciton is no longer visible. 

Figure 1^1 also provides a comparison of cr(ti;) with the 
total density of states: 

7^(u;) = (l/iV)^^(5(c^~i;fc,„) (22) 

k n 

The latter quantity shows apart of the two-domain wall 
continuum, which reaches up to w ~ 0.6V, also the much 
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FIG. 10: Relative optical spectral weights of bound state 
(solid line) and of the first continuum state (dashed) versus 
ti/V calculated for a A'^=26 ring. The binding energy is given 
for small t\ only because of the difficulty to extract this quan- 
tity from a discrete spectrum. The inset shows the kinetic 
energy (solid line) and the result obtained via the sum rule 
(triangles) . 
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larger two particle-hole continuum involving four DW's. 
The latter, however, is not optically active. Comparison 
of D{ijj) and cr(w) indicates a strong enhancement of the 
optical matrix element towards the lower edge of the 2 
DW continuum, i.e., leading to the asymmetric shape of 

cr(w). 

The weight of the lowest (exciton) excitation is com- 
pared in Fig.^lwith the weight of the second lowest exci- 
tation. It can be seen that the decay of exciton weight is 
rather fast with growing ti and parallel to the decrease of 
the exciton binding energy. The inset of Fig. [TUl provides 
a comparison of the kinetic energy due to the quantum 
fluctuations and the integral over the optical conductiv- 
ity as function of ti . The coincidence of the data reflects 
the optical sum rula^S: 

-/ di^a{^) = -Y,5i{-HU. (23) 

In the general case where different hopping processes in 
the kinetic energy i?^j„, I = 1, 2, • • • , reach over differ- 
ent distances 6a the optical sum rule has to be modified 
and is given by an average over the individual kinetic en- 
ergy contributions weighted by the square of the hopping 
distances^^. In all the cases considered thus far, the opti- 
cal sum rule is fully exhausted by the incoherent part of 
(t(u;), thus there is no zero-frequency Drude contribution. 
The increase of total spectral weight with ti naturalUy 
follows from the corresponding increase of kinetic energy. 
The quadratic dependence ~ t\/V at small ti highlights 
that the kinetic energy is due to virtual charge excita- 
tions with energy proportional to V . 

Next we consider the thermal evolution of the opti- 



FIG. 11: Optical conductivity at quarter-filling (p = 0.5): 
(a) t\ — 0.021^ and (b) ti = 0.05V for different temperatures 
T = 0.01,0.10 and 0.20V. 

cal conductivity spectra displayed in Fig. ^2 for the two 
cases ti/V = 0.02 and 0.05. The spectra shown were ob- 
tained for a iV = 26 site cluster and periodic boundary 
conditions. There are two immediately obvious features: 
(i) The disappearance of the prominent exciton absorp- 
tion in the ti = 0.02^ case at high temperature, and (ii) 
the spectral weight transfer from high energy into a low 
energy gapless absorption. Furthermore one also realizes 
a gradual decrease of the total sum rule with increasing 
temperature. 

The appearance of the low energy absorption arises 
from the thermal population of excited states, i.e., ex- 
citon and continuum states, and subsequent transitions 
within the continuum. The evolution of the low-energy 
continuum dictates the temperature dependence of the 
DC conductivity. We shall come back to this point in 
section VI. 



IV. ROLE OF SECOND NEIGHBOR HOPPING 

Edge-sharing chain compounds have the peculiar prop- 
erty that the magnitude of the second neighbor hopping 
matrix element ^2 is larger that the nearest neighbor ma- 
trix element ti . This has the important consequence that 
the Fermi surface topology is changed, i.e., instead of two 
Fermi points there can now be four Fermi points depend- 
ing on the hole or electron concentration. Therefore we 
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FIG. 12: Free fermion dispersion Ek for t2/ti — 2.5 with 
scattering processes Qi = 2kF,i, Q2 = 2fcF,2 and Q3 = 7r/2 
(for p = 1/2) indicated by arrows. 



FIG. 13: Excitation spectrum for = 26, ti = 0.02V and 
two different second neighbor hopping matrix elements t2 = 
0.05 and O.IOV^, showing the strong effect effect of t2 on the 
exciton dispersion with a shift of the minimum to q = n/2. 



have here a qualitatively new situation compared to the 
pure ti case, where the modulations of the CDW arising 
from the Fermi surface instability and of the WL coin- 
cide. At sufficiently large values for t2, i.e., relative to 
the Coulomb interaction strength V, we expect the sys- 
tem to undergo a transition from the Wigner phase into a 
CDW state with different modulation, i.e. now dictated 
not by the classical Coulomb interaction but by 'Fermi 
surface nesting', that is a charge modulation of quantum 
mechanical origin. 

For sufficiently large ^2 there are three relevant scat- 
tering processes Qi = 2kp,i and Q2 = 2kp,2 which are 
in general related to incommensurate modulations deter- 
mined by the ratio t2/ti, and a commensurate modula- 
tion Q3 = 7r/2 as shown in Fig. El Whereas Qi and 
Q2 lead to the opening of gaps at two Fermi points, re- 
spectively, Qs generates gaps at all 4 Fermi points si- 
multaneously. Our numerical results show that the Q3 
modulation is the dominant one, and leads to a further 
doubling of the unit cell. It is remarkable that the mod- 
ulation Q3 coincides with that of the 2A:F-Peierls insta- 
bility at quarter-filling, that is in the usual model with 
nearest-neighbor hopping and spin. Here we shall not 
follow the weak-coupling route further, but investigate 
the transition to the CDW from the strong-coupling WL 
side. 

It is immediately evident that in the perfect charge 
ordered state t2 processes are completely quenched, and 
only through the presence of charge fluctuations intro- 
duced via nearest neighbor hopping ti the t2 processes 
are activated. This is seen in the excitation spectra for 
the Coulomb chain shown in Figs. E| The exact diag- 
onalization results show even for t2 — 0.10 hardly any 
effect on the 2-DW continuum. A small downward shift 
of the 2-DW spectrum is attributed to a broadening of 
the 4DW continuum. 

The exciton, however, is changed in a surprising way. 



it disperses downward, in contrast to the ti case stud- 
ied before. The numerical solution reveals two further 
aspects: (i) the periodicity of the exciton dispersion in- 
dicates nearest neighbor hopping of DW's, quite in con- 
trast to the ti motion where DW's hop over two sites; 
and (ii) the exciton dispersion does not depend on the 
sign of t2. 

For further illustration we present in Fig. 1141 exact di- 
agonalization results for the excitation spectrum of the 
model with an interaction truncated at Imax — 3. The 
domain wall interaction in this case is Ai = 2/3V and 
A,„ — A2 = V for TO > 2. The DW continuum is cen- 
tered around V. The dispersion of the exciton in case 
(a) with ti — 0.02V and t2 = is strongly suppressed 
as compared to the case with long-range Coulomb inter- 
action (cf. Fig. EJ. As we shall see below, this is due 
to the larger splitting A2 — Ai between the exciton and 
the continuum as compared to the model with Coulomb 
interaction 

At finite ^2 = 0.05 (Fig. ^h)) the degeneracy of 
the exciton is lifted and one observes clearly two exci- 
ton branches, one with downward and one with upward 
dispersion. To get some deeper insight in the peculiar 
dispersion of the exciton at finite t2, we turn now to the 
analytical analysis of the excitonic state for the latter 
case, i.e., defined by the DW potential Ai = 2/3V and 
Am = y for TO > 2. 



A. Analytical study of the exciton state 

The analytical calculations presented above can be ex- 
tended to the more realistic case of a finite next-nearest- 
neighbor hopping term t2. For p = 1/2 t2-liopping pro- 
cesses are completely blocked in the ground state with 
perfect charge order. A nearest-neighbor DW pair, how- 
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FIG. 14: Excitation spectrum for a A'^ = 26 site ring at elec- 
tron concentration p = 0.5 and Coulomb interaction trun- 
cated at Imax ~ 3. Numerical results (circles) for ti = 0.02V 
and two different second neighbor hopping matrix elements: 
(a) t2 = and (b) 0.05V. Analytical result for domain-wall 
exciton states shown as solid lines. 



the even and odd sectors in the secular determinant: 



D 





ti{q) UJ2 ti[q) 

h{q) UJi t2{q) 

... q{q) h(q) 







h{q) 0J2 h{q) 
hiq) UJ2 ■•• 

■■. '■. ■■. 



(25) 



The expansion of the secular determinant yields a modi- 
fied equation for the bound states: 



ujip - tiiqY 



|t2(g)|V = 



(26) 



The t2-term leads to a splitting of the degenerate solu- 
tions obtained in the ti-case into a lower and an upper 
branch: 
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FIG. 15: Sketch of 2nd neighbor t2 hopping processes be- 
tween energetically degenerate and symmetrized DW config- 
urations. Such processes occur only in the m — 1 sector 
between symmetric and antisymmetric configurations (due to 
the Pauli principle) and lead to a nearest neighbor motion of 
the center of mass of the DW pair. Signs attached to the right 
brackets indicate even and odd linear combinations. 



ever, can move by one step to the left or to the right 
as a result of a t2-process as indicated in Fig. El For 
the analytical solution it is useful to adopt symmetrized 
configurations as shown in Fig. 1151 

There are only relevant matrix elements ^ t2 between 
even and odd configurations in the m = 1 sector that are 
energetically degenerate: 



'{'4'q,l\Ht2\'4'q,l)+ = t2{q), 



(24) 



Ee^,l{q) ^ Ai-|<2(g)| 
Ee^^uiq) - Ai + \t2iq)\ 



A2-Ai + |i2(9)r 
A2 - Ai - \t2{q)\' 



(27) 
(28) 



where t2{q) — 2it2sin{q). This matrix element connects 



The coupling t2 shifts the exciton minimum in the lower 
branch to 7r/2. The exciton dispersion does not depend 
on the relative sign of ^2- At a threshold value t2,cr 
q — tt/2 excitons condense and lead to the new state 
controlled by the kinetic energy. The upper branch is a 
physical solution only as long as it does not touch the 
DW continuum. 

In Fig. 1141 we provide a comparison of the analytical 
solutions for the exciton dispersions, Eas. H27l28l ). with 
exact diagonalization data for a = 26 site cluster in 
the case of Coulomb interaction truncated at Imax — 3. 
As can be seen from Fig. 0]in this case Ai — 2/3F and 
Am = A2 = y for m > 2. The analytical exciton dis- 
persions calculated from Eas. H27l28l ) provide a good de- 
scription of the numerical values. The small deviations 
are mainly due to a slight downward shift of the two-DW 
continuum in the numerical calculation, that results from 
its interaction with the 4-DW continuum which is much 
broader in the t2 = 0.05 case. This also implies a small 
downward shift of the DW-exciton. 

It is also instructive to compare Fig. I14r b) calculated 
for the truncated interaction with Fig. Il^f a) which was 
obtained using the full Coulomb interaction. In the latter 
case only the lower exciton branch can be seen and its 
dispersion is slightly larger than in Fig. I14r b) although 
the spectra have been determined for the same hopping 
parameters. These differences result from the different 
attraction of domain walls A2 — Ai in the two cases, 
which enter in the 3rd term on the r.h.s. of exciton dis- 
persion in Ea. (|27l ). 



B. Structure factor and exciton instability 

In the previous discussion of the electronic structure 
we have seen that the exciton state will get soft at about 
<2 = O.lSy and hence the WL state should get unstable 
and a new ground state with momentum g ^ 7r/2 should 
appear. In the following we analyse the change of the 
static charge structure factor. Figure ^1 shows the de- 
pendence of the charge structure factor N{q), Ea. ((T^ . 
on the size of t2 for fixed ti = 0.021^. The results are 
obtained for a iV = 10 site cluster where we included the 
full Hilbert space in order to show the complete break- 
down of Wigner order at large t2 ■ 
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FIG. 17: Excitation spectrum for L — 26, ti — 0.02V and two 
different second neighbor hopping matrix elements t2 = 0.15 
and 0.20 V indicating the exciton instability of the Wigner 
lattice. While for t2 = 0.15 the degenerate ground states are 
at q = 0, n, the ground state momenta for t2 — 0.2 are near 
q = 7r/2. 



We note, that our program was developed for iV = 47i + 2 
site rings {n — 1, 2, • • • ) that do not have 7r/2 as allowed 
momentum, therefore there are two degenerate ground 
states at (7 = 7r/2 ± 2'k/N . 

The critical value t2,cr can also be estimated from 
Ea. (|77|) by setting Eex,i{q) = at g = 7r/2 and using 
the parameters for the DW interaction appropriate for 
the model with Coulomb interaction. 



FIG. 16: Charge structure factor N{q) as obtained by ED for 
a 10 site ring with Coulomb interaction {V = 1) and p = 0.5 
for U = 0.02 and different values for t2 (T = O.OIV). Straight 
Unes are guides to the eyes; the -I- symbols in the data for 
t2 = 0.16 also indicate the allowed q points. Wigner order is 
unstable beyond the critical value t2^ — 0.155. 

It is remarkable that up to the value t2 < 0.14 the 
ground state correlations remain unchanged with N{q) 
peaked at momentum tt. This reflects the blocking of 
t2 hopping processes in the state with alternating charge 
order. Then in the narrow range 0.15 > t2 > 0.16 there 
is a sudden change of N{q) indicating a level crossing. 
Beyond t2 ~ 0.16 the charge correlations are determined 
by the new state with the maximum of N{q) near tt/2. 
The broad shape of N{q) is reminiscent of that of a ID 
Fermi gas, i.e., indicating that the CDW modulations 
in this phase are rather weak. It has been checked by 
a calculation of the charge correlation functions in real 
space, that there are still significant charge modulations 
in the system consistent with a. q = tt/2 CDW. 

The excitation spectrum corresponding to the transi- 
tion from the WL to the 7r/2 CDW is shown in Fig. ITTfa). 
It clearly reveals the quasi-degeneracy of the WL ground 
states aX q — 0{tt) and the soft exciton at g = 7r/2. At 
larger values for t2 the state emerging from the exciton is 
the new ground state of the system, as seen in Fig. llTt b). 



C. Optical conductivity 

The effect of second neighbor hopping t2 on the opti- 
cal conductivity a{uj) is displayed in Figs. irHT a.b) for the 
model with Coulomb interaction. As a consequence of 
the blocking of t2 processes in the state with alternating 
Wigner charge order at p — 1/2 the optical spectra are 
only slightly modified at low temperature, and almost 
coincide with those obtained for the model with only 
nearest-neighbor hopping Fig. 1111 A remarkable change, 
however, is the disappearance of the q = exciton, which 
was very pronounced in the spectrum for ti = 0.02 and 
t2 = in Fig. lllf a') at low temperature. This feature has 
disappeared after switching on t2 = 0.05 in Fig. HHf a). 
For these parameters there is no bound state at q — 0(tt) 
as can also be seen from the corresponding energy level 
diagram in Fig. Iisr a) . 

As an optical experiment involves only vertical transi- 
tions only charge excitations at q — and tt are probed at 
low temperature, i.e., the downward dispersing exciton is 
invisible. This changes, however, when the exciton states 
get populated by thermal excitations. As a consequence 
there are marked changes in the spectra at higher tem- 
perature, which can be traced back to the different dis- 
persion of the exciton state at finite t2 ■ The pronounced 
structure near uj = 0.2 in Fig. llSt a) at T = O.IT^, which 
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FIG. 18: Optical conductivity at quarter-filling (p = 0.5) for 
different temperatures: (a) ti = 0.02, t2 — 0.05 (Lorentzian 
broadening 7 = 0.05V), and (b) h = 0.02, t2 = 0.10 . The 
in-gap ansorption appearing near ui ~ 0.19V at finite temper- 
atures in (a) is attributed to a transition from the thermally 
excited exciton in the vicinity of g = tt to the domain wall 
continuum. A similar feature appears in (b) at uj ^ 0.27V. 



U — V2 with intensity ~ But also in the case ^2 *C ti 
transitions into the upper Hubbard band are expected 
due to charge fluctuations resulting from ti processes. 

V. DC-CONDUCTIVITY 

A central experimental quantity to compare our finite 
temperature results with is certainly the DC-conductivity 
cr(T). In Fig. which displays the temperature de- 
pendence of cr^Lu), we observe the emergence of a low- 
frequency continuum whose intensity grows with increas- 
ing temperature. These changes are accompanied by a 
spectral weight transfer from high to low energy. The 
low-energy excitations arise from transitions within the 
2 DW continuum. We note that at higher temperatures, 
i.e., near the melting temperature of the WL and above, 
also 4 DW excitations do contribute substantially to the 
low-frequency absorption. This yields a very dense low- 
energy spectrum even for small systems. It is certainly 
suggestive that this low-frequency continuum contains 
the information about the DC-conductivity. Yet as we 
are dealing with finite systems the case is not that sim- 
ple as the zero frequency limit cannot be determined in 
a straightforward way. In fact the analysis of the low- 
frequency part of cr(w) reveals a pseudogap ^ 8t/N which 
scales inversely with the size of the system N. 

Similar pseudogap behavior and finite size effects in 
(j{w) have also been observed recently by Prelovsek et 
alJ^ in a study of the ID t-V model. In their careful 
study Prelovsek et al. arrived at the conclusion that after 
finite size scaling the results for frequency dependence 
could be compatible with a normal and featureless shape 
of (T(tj) as found by a frequency moment analysis. 



is already seen at T = 0.05V as a weak in-gap excitation, 
stems from transitions between thermally excited exci- 
tons near q — tt/2 and final states at the upper edge of 
the 2DW continuum ( cf. Fig.^Ja)). This conclusion is 
based on a careful study of the size of the corresponding 
matrix elements in the current-current correlation func- 
tion. In the spectra of Fig. llSf b) for a twice as large 
t2 value this absorption has shifted to higher energies 
(lu ~ 0.27V) and appears now as a small structure on 
the shoulder of the main peak, which is due to transi- 
tions into the 2DW continuum at g = 0. 

Thus the observation of an in-gap absorption in the 
optical conductivity at finite temperature can provide 
valuable information about the position of the minimum 
of the exciton dispersion in the middle of the Brillouin 
zone, which may allow to determine the value of t2 from 
experiment. 

Finally we recall that the higher transitions into the 
upper Hubbard band with energy ~ U are not contained 
in the spinless fermion model we study here. Yet they 
are contained in the Hubbard- Wigner model. For t2 > 
ti and p = 1/2 these transitons are expected at energy 



Thus we proceed as follows: We assume that in the 
thermodynamic limit the level spacing vanishes and and 
the low frequency part of the spectrum can be expressed 
by a Drude form for the real-part of the conductivity 



a{uj) 



1 + {uOTf 



(29) 



where r w 1/ii is determined by the energy scale of the 
domain-wall continuum. The spectral weight of the low- 
frequency part of cr(a') as determined by exact diagonal- 
ization 



/(wo) 



a{uj)<kij, 



(30) 



with Wo = 4ti is then used to determine the DC- 
conductivity (j{T) with help of Ea. l(^ . 

Results for the DC-conductivity obtained from a 26- 
site ring at doping p = 0.5 are presented in Fig ^] for two 
different hopping matrixelements. The numerical data 
for ti — 0.02 and 0.05V (t2 = 0) show similar behavior. 
The DC conductivity reveals activated behavior below 
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FIG. 19: (a) Temperature dependence of the DC- 
conductivity for p = 0.5 determined from the integrated low- 
frequency spectral weight for ti — 0.02 and 0.05 {t2 = 0). (b) 
Arrhenius plot revealing the activated behavior at low tem- 
perature. An activated behavior for the carrier density in 
combination with a temperature dependent mobility ~ T~°, 
Ea. H^l^ . provides a rather good fit of the numerical data over 
the whole temperature range. 



and basically T-independent behavior above the mehing 
temperature of the WL. The conductivity for ti = 
0.02V is smaller, reflecting the proportionality of the 
current-current correlation function. 

At small hopping ti = 0.02^ the numerical data for 
cr(T) is well described over 5 decades by a single acti- 
vation energy dependence a(T) ~ cexp—Ea/kT, with 
Ea — 0.30V. However, the sometimes used relationSSiZi 

a(T) ~ cT-" exp {-Ea/kT), (31) 

which involves in addition the temperature dependence 
of the mobility fi{T) ~ T~", improves the fit when ap- 
proaching the saturation (melting) regime at high tem- 
perature as shown in Fig^jby dashed lines. Using this 
relation we obtained for ti — 0.02 (0.05)1^ the activa- 
tion energy Ea ^ 0.39(0.32)F, respectively, and a ~ 
2.4(1.95). Thus the mobility /i(T) of the carriers (domain 
walls) decreases strongly with increasing temperature, as 



one may have expected. The activation energies Ea de- 
termined from a{T) are consistent with the energy gaps 
in the corresponding energy level diagrams in FigsEIa,b). 
Yet the activation energies Ea determined from a{T) are 
larger than expected from the relation Ea — Eg/2 which 
applies for usual semiconductors. Whether this discrep- 
ancy originates from the fact that the charge carriers in 
the WL are domain walls with fractional charge and not 
usual electron- and hole-like quasiparticles remains un- 
clear and deserves further study. 

While the relation Eq. (|31|l provides a quite satisfactory 
description of the numerical data over the full tempera- 
ture range, it is nevertheless far from perfect. As one can 
see, the fit curve in Fig tTW b) lies below (above) the nu- 
merical data for temperatures above (below) the melting 
temperature T^, respectively. This deviation possibly 
reflects the strong change of the structure factor near T^, 
(see inset of FiglSj). Here our aim was to keep the num- 
ber of parameters as small as possible and therefore we 
have not tried to add such a Ea{T) correction term that 
would involve Tm and further parameters. 

Recent measurements of the DC-conductivity of 
Na3Cu204 and NasCusOio compounds show the same 
trends: (i) an Arrhenius behaviour below , and (ii) a 
saturation of DC conductivity above the melting tran- 
sition Tjn ~ 455(540)K, respectiveljiSiZ^. Moreover, the 
experimental conductivities show only a small disconti- 
nuity at the melting temperature T„i. Although the ab- 
sence of a discontinuity in the theoretical curve can be 
attributed to the finite system, which is not expected 
to display a phase transition, it is nevertheless remark- 
able that also the experimental data does only show a 
weak discontinuity at the melting temperature of the 
WL. The measured DC-conductivity of the compound 
Sri4Cu2404i was also found to be described by an Ar- 
rhenius law with an activation energy Ea ~ 0.12 eV. 
More recent studies found a crossover between two expo- 
nential regimes, with Ea ~ 0.12 eV for the low T regime 
up to about 170 K but different values 0.18 eVi&, 0.22 
eV^s and 0.27 eV^i for the high-T regime. Up to 400 
K no saturation was observed for this compound. As 
the transport at high temperatures is probably due to 
the chains and the ladders in these compounds, a direct 
comparison with our results is ruled out. 



VI. CONCLUSIONS 

In summary, we have investigated the charge excita- 
tions of a ID generalized Wigner lattice, expected to be 
realized in edge-sharing Cu-0 chain systems and also 
in some organic chain compounds, starting from the 
Hubbard- Wigner model with long-range Coulomb inter- 
actions V/l among electrons. A central aim was to gain 
insight into the spectral structure of the optical conduc- 
tivity a{uj) and its temperature dependence. 

We have found that: (i) a{u!) is determined by a highly 
asymmetric spectrum due to two domain- wall excitations 
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with an energy gap Eg ~ V/2 — Ati and a width of ~ 
8ti. The asymmetric form of this spectrum may serve 
as another fingerprint to detect the WL and may also be 
employed to determine the parameters of the model from 
experiment. 

(ii) For ti <ti V excitons with a dispersion given by 
Ea. (|18|l form and show up as strong absorption peak 
in the optical conductivity. The excitons are due to an 
attractive potential between domain- wall pairs ^ V/81, 
whose prefactor reflects the fractional charge of the 
DW's, and result as a consequence of the long-range (re- 
pulsive) Coulomb interaction between electrons. The ap- 
pearance of excitonic states, which arise from the effec- 
tive attraction between the fractionally charged DW's, 
was to our knowledge not noted before. In contrast in 
the frequently studied model with only nearest neighbor 
interactions there is no bound state because of the ab- 
sence of attractive interactions between DW pairs. How- 
ever, in that case interchain interactions may provide an 
alternative mechanism for confinement, as noted recently 
in work by Bhaseen and Tsvelik--. 

(iii) Edge-sharing chain compounds have the unusual 
property that the magnitude of the second neighbor hop- 
ping matrix element t2 is larger than the nearest neighbor 
matrix element ti. While t2 hopping processes are frus- 
trated in the classical WL state with alternating charge 
order, these processes surprisingly contribute strongly to 
the exciton dispersion. They lift the degeneracy of the 
exciton state and the lower branch, c.f. Ea. H27|l . leads to 
an exciton instability at about t2,c ~ 0.16^. The CDW 
state beyond t2,c has a modulation period twice as large 
as that of the WL. The charge modulation is weak in this 
CDW state as infered from the calculation of the static 
charge structure factor. 

(iv) Interestingly the optical conductivity at finite tem- 
perature reveals an in-gap absorption which reflects the 
transitions between the soft exciton near q = 'k/2 and 
the domain wall continuum. Thus this in-gap absorption 
may provide a further identification of the WL state and 



also allow for an independent experimental determina- 
tion of the matrix element ^2- 

(v) Moreover, we have calculated the temperature 
dependence of the DC-conductivity of the generalized 
Wigner lattice from the low-frequency absorption. The 
data for a{T) shows a crossover from activated behav- 
ior at low temperature to a basically temperature inde- 
pendent conductivity at high temperatures. It turned 
out that a{T) can be described over the full temperature 
range by an activated behavior characterized by an acti- 
vation energy Ea and a temperature dependent mobility 
^ 2^-". This implies a strong decrease of the mobility of 
domain walls with increasing temperature. 

Finally, doped edge-sharing compounds provide a 
unique opportunity to study the competition between 
two entirely different states, the classical WL dictated 
by the long-range Coulomb interaction and the CDW of 
quantum mechanical origin, i.e., resulting from a Fermi 
surface instability. These materials highlight the impor- 
tance of long-range Coulomb interaction in strongly cor- 
related systems, — and provide a one-dimensional test 
ground for the study of charge stripe formation. We 
have analysed here the charge dynamics and aspects of 
transport for a generalized ID Wigner lattice in its most 
simple realization, namely at quarter-filling [p — 1/2). 
Other commensurabilities are more complex and show 
a hierachy of different charge excitations. Work along 
these lines is in progress, as well as work on the effect 
of spin-charge coupling particularly on the low-energy 
charge response. 
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